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Abstract 

We consider spherically symmetric black holes in generic Lovelock gravity. Us- 
ing geometrodynamical variables we do a complete Hamiltonian analysis, in- 
cluding derivation of the super-Hamiltonian and super-momentum constraints 
and verification of suitable boundary conditions for asymptotically flat black 
holes. Our analysis leads to a remarkably simple fully reduced Hamiltonian 
for the vacuum gravitational sector that provides the starting point for the 
quantization of Lovelock block holes. Finally, we derive the completely reduced 
equations of motion for the collapse of a spherically symmetric charged, self- 
gravitating complex scalar field in generalized flat slice (Painleve-Gullstrand) 
coordinates. 
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1 Introduction 



The origin of the four-dimensionality of the present universe is one of the most funda- 
mental problems in gravitational physics. One possible explanation is that only the four- 
dimensional universe is stable in some physical sense and therefore chosen at the moment 
of its creation. Another possibility is that after the creation of a higher-dimensional uni- 
verse, the extra spatial dimensions other than our perceived three-dimensional space are 
compactified by some mechanism. Einstein's general theory of relativity allows us to study 
this fundamental problem by setting the number, n, of spacetime dimensions as a tunable 
parameter. It is therefore possible to look for critical values of n beyond which spacetime 
properties change drastically. Such analyses can give us valuable insights about the origin 
of our four dimensional universe. 

In this context, it is reasonable to focus on black holes because they are fundamental 
objects that encode many key features of the gravitational interaction. In four dimensions, 
asymptotically flat stationary black holes are characterized by a small number of parameters 
such as mass or angular momentum, somewhat analogous to atomic properties in chemistry. 
This is a consequence of the black-hole uniqueness theorem asserting that the Kerr- Newman 
black hole is the unique asymptotically flat stationary and rotating black hole with a 
connected horizon in the Einstein-Maxwell system. (See [1] for review.) It is important to 
note that this uniqueness theorem is not valid in higher dimensions j2], as explicitly shown 
in five dimensions, for example, by the existence of two distinct asymptotically flat black 
objects with the same mass and angular momentum but with different horizon topology [31 
H]. This fact already demonstrates one special feature of four-dimensional spacetime. n = 4 
is further singled out as a critical value in the framework of general relativity because 
asymptotically flat vacuum black holes do not exist in n < 4 dimensions [5]. 

When studying higher dimensions, it must be remembered that general relativity is not 
the only natural extension of four- dimensional Einstein gravity. General relativity is a 
quasi-linear second-order theory, which ensures the well-definedness of the initial value 
problem and the absence of ghosts. In 1971 Lovelock showed that in higher dimensions 
general relativity is just a special case of the most general class of theories satisfying these 
property. These more general theories are collectively called Lovelock gravity [6]. The 
Lovelock Lagrangian consists of a sum of the dimensionally extended Euler densities in 
which the cosmological constant and the Einstein-Hilbert terms appear as the zeroth- and 
the first-order terms, respectively. Just as the Einstein tensor trivially vanishes in two- 
dimensional spacetime, the second-order Lovelock Lagrangian (called the Gauss-Bonnet 
term) becomes purely topological in four spacetime dimensional and does not contribute 
to the field equations [TJ. As a consequence, Lovelock gravity reduces to general relativity 
with a cosmological constant in four dimensions. 

Motivation for studying Lovelock gravity is also provided by the fact that string/M- 
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theory [8] requires the existence of extra spatial dimensions whereas quantum theory sug- 
gests the need to add higher-curvature terms. Moreover, string theoretic arguments [H] 
suggest that quadratic Lovelock gravity appears in the low-energy limit for strings propa- 
gating in curved spacetime. 

For the above reasons, Lovelock gravity has been extensively investigated, with emphasis on 
the similarities to and difference from general relativity. (See [10] for review.) However, in 
comparison with its classical aspects, the quantum theory is poorly understood at present, 
despite the fact that deep issues such as the origin of the four-dimensional universe can 
only be answered in the quantum context. In four-dimensional general relativity, Kuchaf 
presented an elegant geometrical framework for studying the physical phase space and 
quantization of spherically symmetric vacuum black holes [11] . Using Kuchaf geometrody- 
namics as a foundation, Louko and Makela were able present a rigorous quantization of the 
Schwarzschild black hole spacetime, including a construction of all self-adjoint extensions of 
the Hamiltonian and derivation of the semi-classical area spectrum [12]. They found that 
the area/entropy spectrum was equally spaced in the semi-classical limit, in agreement with 
early speculations of Bekenstein and Mukhanov based on the thermodynamic properties 
of black holes [13]. An equally spaced area spectrum was also obtained for Schwarzschild 
black holes using a variety of different techniques [14J . This result is perhaps not surprising 
because there is only one length scale in the system, namely the Planck length. In con- 
trast, there is more than one length scale in quantum Lovelock gravity because the couping 
constants to each order of the Lovelock Lagrangian are dimensionful. As a consequence, 
the entropy of Lovelock black holes is no longer equal to 1 /4 the area [153 • ^ * s therefore 
of great interest to see what quantum spectrum emerges for the both the area and entropy. 

It is important to mention two key features of general relativity that were crucial to Kuchaf 's 
analysis. First the existence of Birkhoff's theorem implies that the reduced phase space 
is finite dimensional (two dimensional in the case of Schwarzschild black holes). Secondly, 
Kuchaf identified the Misner- Sharp mass [16] in spherically symmetric vacuum spacetime 
and its conjugate momentum, the Schwarzschild time separation at infinity, as the physical 
phase space variables. In more general situations (i.e. non- vacuum), the Misner-Sharp mass 
is known as the best quasi-local mass in the spherically symmetric spacetime. It satisfies 
the requisite monotonic and positivity properties and converges to the Arnowitt-Deser- 
Misner (ADM) mass at spacelike infinity in asymptotically flat spacetime [IT]. Fortunately, 
Lovelock gravity also possesses both these key features. Bikrhoff's theorem in Lovelock 
gravity asserts that the spherically symmetric vacuum solution is uniquely determined 
under certain conditions [18]. The corresponding Schwarzschild- Tangherlini-type vacuum 
solution was obtained by Zegers [18J. In addition, a natural counterpart to the Misner- 
Sharp mass has been defined in Lovelock gravity [19, 20J. 

The purpose of the present paper is to provide a framework to study quantum aspects of 
spherically symmetric black holes in Lovelock gravity. We provide a comprehensive analysis 
that goes far beyond the initial presentation of our results in [21]. In particular, we use the 



4 



geometrodynamical formulation of Kuchaf to do a complete Hamiltonian analysis, including 
derivation of the super-Hamiltonian and super-momentum constraints and verification of 
suitable boundary conditions for asymptotically flat black holes. Our analysis leads to a 
fully reduced Hamiltonian that is just as simple as that of Kuchaf. As a specific application, 
we also derive the fully reduced equations of motion in flat slice coordinates for the collapse 
of a charged scalar field, including Lovelock gravitational as well as electromagnetic self- 
interactions. 

We note that the Hamiltonian analysis for full Lovelock gravity was first considered by Teit- 
elboim and Zanelli [22]. Their result was rather formal in that an explicit parametrization 
of the phase space is not provided. (See also [231 121] •) For the case of spherical symme- 
try, the geometrodynamics [11] of five- dimensional Einstein-Gauss-Bonnet (i.e. quadratic 
Lovelock) gravity was worked out by Louko et al [22]) while the Hamiltonian analysis of 
higher-dimensional Gauss-Bonnet gravity coupled to matter was recently done in [26] . Our 
analysis is done for generic Lovelock gravity in arbitrary dimensions. 

In the following section, we present our system, including action and spherically symmetric 
solutions. Section III derives a dimensionally reduced equivalent two-dimensional action 
that is the starting point of our analysis. It also reviews the geometrodynamics of Kuchaf 
in a general dynamical setting. In Section IV, we perform the Hamiltonian analysis for 
general relativity, in terms of both the standard ADM and geometrodynamical variables. 
The generalization to Lovelock gravity is presented in Section V. The contributions of 
matter fields are discussed in Section VI, while concluding remarks and discussions appear 
in Section VII. 

Our basic notation follows [27]. The convention for the Riemann curvature tensor is 
[V p , V ff ]y = W i vpa V v and 1Z^ U = 1Z P fipu . The Minkowski metric is taken as diag(— , +, • • • , +), 
and Greek indices run over all spacetime indices. We adopt the units in which only the 
n-dimensional gravitational constant G n is retained. 



2 Preliminaries 



2.1 Symmetric spacetimes in Lovelock gravity 

The action of the gravitational system is written as 

I — Im + IdM, (2-1) 

where Im is the dynamical term and Iqm is the boundary term. In general relativity, 
Im is the Einstein-Hilbert action and Iqm is the Gibbons-Hawking- York boundary term. 
(See [2H] for the boundary term in general Lovelock gravity.) In the present paper, we 
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consider Lovelock gravity in n(> 4)-dimensional spacetime, of which dynamical term in 
the action is given by 




where K n := y/8nG n . Our notation basically follows [19]. «( p ) is the coupling constant for 
the pth-order Lovelock Lagrangian with dimension (length) 2 ^ -1 ) and we assume /t 2 > 
without any loss of generality. The 5 symbol denotes a totally anti-symmetric product of 
Kronecker deltas, normalized to take values and ±1, defined by 

••-<%• (2-4) 



The gravitational equation following from this action is given by 

Q \lu — f^nTf^u, (2-5) 

where v is the energy-momentum tensor for matter fields obtained from i ma tter and 

[n/2] 

^:=5><p)Gjg, (2-6) 

p=0 

nvip) — \ xmi-VpCy-Cpn P iai ...i? p ^ (9 7\ 

The tensor G$ obtained from C( p ) contains up to the second derivatives of the metric and 
G$ =0is satisfied for p > [(n + l)/2]. 

In the present paper, we consider the n(> 4)-dimensional warped product spacetime 
(M n ,g^ u ) ~ (M 2 ,g AB ) x (K n ~ 2 ^ ab ) with the general metric 

g^{x)dx»dx v = g AB {y)dy A dy B + R{y) 2 lab {z)dz a dz\ (2.8) 

where g AB is an arbitrary Lorentz metric on (M 2 ,g AB ) and R(y) is a scalar function 
on (M 2 , g AB ). lab is the metric on the (n — 2)-dimensional maximally symmetric space 
(K n ~ 2 ,'y ab ) with its sectional curvature k = 1,0,-1. We note that the results in the 
present paper are valid for k = with p = by setting k p = 1. We introduce the covari- 
ant derivatives on spacetime (M n 1 g iiu ) ) the subspacetime (M 2 ,g AB ) and the maximally 
symmetric space (K n ~ 2 ,^ ab ) with 

V p ^ = 0, D F g AB = 0, D f g ab = 0. (2.9) 
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The most general energy-momentum tensor T^ v compatible with this spacetime symmetry 
governed by Lovelock equations is given by 

T, u dx»dx» = T AB (y)dy A dy B + p(y)R 2 lab dz a dz b , (2.10) 

where T AB (y) and p(y) are a symmetric two-tensor and a scalar on (M 2 , g AB ), respectively. 

The generalized Misner-Sharp mass in Lovelock gravity is defined by 

(n - <2W {k) [n/2] 

M := [n -f^g * {p) R n - l - 2p [k - (DR) 2 f, (2-11) 

n p=0 



(n — 3)\a( v ) 

S « := (n-ll^l ' (2 - 12) 



where (Di?) 2 := {D A R)(D A R) [31]. The constant represents the volume of (K n ~ 2 , j ab ) 
if it is compact and otherwise arbitrary positive. M reduces to the ADM (Arnowitt-Deser- 
Misner) mass at spacelike infinity in the asymptotically flat spacetime. In terms of M, some 
components of the Lovelock equation are written in the following simple form [T9l [29] : 

D A M =V^%R n - 2 (t a b (D b R) - T B B (D A R)j . (2.13) 



2.2 Vacuum solutions 



In the vacuum case {T^ v = 0), Eq. (12.131) shows that M is constant. The maximally 
symmetric solution, namely Minkowski, de Sitter (dS) or anti-de Sitter (AdS) solution, gives 
M = 0. The maximally symmetric spacetime may be given in the following coordinates: 

ds 2 = —{k — \r 2 )dt 2 + ^- + r 2 lab dz a dz\ (2. 14) 

k — Xr 2 

where A := 2A/[(n — l)(n — 2)] and A is the effective cosmological constant, which is 
determined by the following algebraic equation: 

[n/2] 

= ^a (p) A p =: v(X). (2.15) 

The Minkowski vacuum (A = 0) is possible only if a™ = 0. Since Eq. (" 12 . 1 5[) is a higher-order 
polynomial, there can be multiple values of A. We call the vacuum A = Ai non- degenerate if 
(dv/d\)(Xi) 7^ holds. A simply degenerate vacuum is characterized by (dv/dX)(Xi) = 0, 
while a doubly degenerate vacuum is characterized by (dv/dX)(X\) = (d 2 v/dX 2 )(Xi) = 0. 
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In a similar manner, a gth-order degenerate vacuum is defined by (d s v/d\ s )(\i) = for 
s = 1,2, ■ ■ ■ ,q, where q < [(n — 3)/2] is satisfied because of <S( n /2) = for even n. 

The Schwarzschild-Tangherlini-type vacuum solution in Lovelock gravity [H] is given by 

dr^ 

ds 2 = -f(r)dt 2 + — + r 2 lab dz a dz b , (2.16) 
f(r) 

where the metric function /(r) is determined algebraically by 

[n/2] 

M = J2^( P )r n - 1 - 2p (k-f(r)Y. (2.17) 

p=0 

M is related to the constant generalized Misner-Sharp mass as M := 2K 2 l M/[(n — 2)V^ 2 ]- 
This class of vacuum solutions in higher-order Lovelock gravity was first obtained by Boul- 
ware and Deser [30] in the quadratic theory. (See [3T] for further discussions.) The above so- 
lution reduces to the ones found in [321 [33] in the case where the coupling constants are cho- 
sen such that the theory admits a fully degenerate maximally symmetric vacuum. Birkhoff 's 
theorem in Lovelock gravity asserts that, in the case where (DR) 2 := (DaR)(D a R) ^ 
and the spacetime is of the C 2 -class, there is a unique vacuum solution as long as the theory 
does not admit degenerate vacua [18]. (See also [19] for more general case.) One of the 
purposes of the present paper is to derive the formulae to quantize a Lovelock black hole 
described by the above solution. 

If the theory admits degenerate vacua, there may be more vacuum solutions. If the theory 
admits simply degenerate vacua, the following is also a vacuum solution: 

ds 2 =-{k- \r 2 )e 25 ^ r) dt 2 + - + r 2 lab dz a dz b , (2.18) 

k — Ar 2 

where S(t,r) is an arbitrary function and A takes the value for the degenerate vacuum. 
(This solution was first obtained properly in quadratic Lovelock gravity by Charmousis 
and Dufaux [M]-) If the theory admits doubly degenerate vacua, there is another vacuum 
solution where the two-dimensional portion (M 2 ,qab) is totally arbitrary [19] . If one re- 
moves the assumption of C 2 -differentiability of the spacetime, then more vacuum solutions 
exist (33 Eg]. 
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3 Dimensionally reduced action 



3.1 Covariant form 



In the symmetric spacetime under consideration, the system may be described by the 
effective two-dimensional action: 



(3.1) 
(3.2) 



1(2) = hi + I&M- 

The dynamical term Im is written as 

Im = J dy°L[g AB ,R] = J df J dfC^AB^R], 

where y° is a timelike coordinate on (M 2 , <7ab)- Here the Lagrangian L and the Lagrangian 
density C are functionals of the metric functions, which are determined up to a total 
derivative. The main purpose of geometrodynamics is to find canonical variables (that are 
functionals of the metric functions) to provide a tractable form and transparent physical 
meaning for L. 

Our first task is to derive a tractable tensorial form of Im- For symmetric spacetimes under 
consideration, the action reduces to 



n J p= o 



(3.3) 



where gp) '■— det(gAB) and the dimensionally reduced pth-order Lovelock term £( p ) is given 
from expressions (2.19) and (2.20) of [19J as 

,fk-(DR) 2 Y , n 2 R fk-(DR) 2 ^^ 1 

(n - 2p)(n - 2p - 1) K — -L - 2p(n - 2p)- 1 



C 



(p) 



(n-2)! 
(n- 2p)\ _ 



R 2 



R 



R 2 



+ 2p{p _ 1) {D*R?-(D a DbR){D b DaR) ( k-(DR) 2 y- 2 | p gj (k - (DR) 2 ^ 



R 2 



R 2 



R 2 



(3.4) 

where is the Ricci scalar on (M 2 ,gAB) and D 2 R := D a DaR- At a glance, there is a 
non-minimal coupling between (DR) 2 and ^TZ in Such a Lagrangian is not tractable 
to perform the canonical analysis. As proven in Appendix I A. 1\ we can write it, up to total 
divergences, without such a coupling: 



C 



(p) 



(n-2)! 
(n-2p)\ 



.(2) 



pk p - l nR 2 ~ 2p + pR 



2-n 



D A {R n - 2p )D A {{DRf 
{DR) 2 



kT 1 — (k— (DR) 2 )*- 1 



+ (n - 2p)(n - 2p - l)!^(k - (DR) 2 ) P + 2pk p -\DR) 2 ^R~ 2p 
This is a key result of our paper and the starting point of our canonical analysis. 



(3.5) 
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3.2 ADM form 



We are going to write down the action ( 13. 5ft by adopting the following ADM coordinates 
(t,x) on (M 2 ,g AB ): 



ds\ 2) = g AB dy A dy B = -N(t, xfdt 2 + A(i, x) 2 (dx + N r (t, x)dt) 2 . 



(3.6) 



Now canonical variables are N, N r , A, and R and their momentum conjugates are respec- 
tively written as P/v, P/v r , Pa, and Pr. In the present paper, a dot and a prime denote a 
partial derivative with respect to t and x, respectively. The metric and its inverse are 



g tt = -(N 2 -A 2 N 2 ), g tx =A 2 N r , g x 



A 5 



g u = -N- z , g tx = N r N~ 
while y/—g(2) is given by 



N~ 2 A- 2 {N 2 - A 2 N 2 ) 



y/-9m = NA. 
For the later use, we compute the following quantities: 

F : = (DR) 2 

= -y 2 + A- 2 R'\ 
y^g~^D 2 R = - d t (Ay) + d x (AN r y + A^iVP') 



where y is defined by 



y := N-\R-N r R'). 



(3.7) 
(3.8) 



(3.9) 



(3.10) 
(3.11) 

(3.12) 



We also need the following relationship: 

D A {R^)D A {{DRf) = (n-p)R n ^- 1 ^yF+(^- 2 +^y)R'). (3.13) 
Using this result, the action ( 13. 5 \} is written in the following simple form: 



M = — 2< — ^ / y v-9w; 

— p(n — 2p) 



a 



(p) 



p=0 " 
jjm— 2p— 1 



n — 2p) 



(2) 

pk p - l KR n - 2p 



y x {hT 1 - (k - F)^}|a^ - (AN r y + A-'NR')^ 



+ {n-2p){n-2p-l)\^k- Ff + 2pk p - l F^R 



n-2-2p 



(3.14) 
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The first term is the two-dimensional gravity non-minimally coupled scalar field R, which 
is essentially the same as the general relativistic case. This term can be explicitly written 
down in terms of the canonical variables using 

y^^R^n = - 2AT 1 [[R n - 2p )'N r - d t (R n - 2p ^j {N' r A + N r A' - A) 

- 2N^(R n - 2p )"A- 1 + (iT-*) / (A- 1 )^ +d t {- ■■)+d x (---). (3.15) 

Based on the action (13.141) . we will perform the canonical analysis in the subsequent sec- 
tions using geometrodynamical phase space variables. We therefore now review briefly the 
geometrodynamics of Kuchaf [11] . 



3.3 Geometrodynamics 

The metric (I3.6P may be written in the generalized Schwarzschild form in terms of the areal 
coordinates as 

ds 2 {2) = —F(R, T)e 2a{R ' T) dT 2 + F(R, T^dR 2 . (3.16) 
The generalized Misner-Sharp mass M is then given by 

(r> -2)V {k) [n/2] f \ p 

M{R, T) = \ - 2 V a {p) R n - l ~ 2p ( k — F(R, T) ) . (3.17) 

This implicitly gives the functional form F = F(R,M). However, there is no one-to-one 
correspondence between F and M unless all the coupling constants are non-negative. 

To see the relation to the ADM form (13. 6 j) we use the coordinate transformations T = 
T(t,x) and R = R(t,x), to write the metric (I3.16P as 

ds 2 2) = - {Ff 2 e 2a - F- x R 2 )dt 2 + 2{-FtT'e 2a + F^RR^dtdx 

+ (-FT' 2 e 2a + F- l R' 2 )dx 2 . (3.18) 

Comparing with the ADM form, we identify 

FT 2 e 2cr - F- X R 2 =N 2 - A 2 N 2 , (3. 19) 

—FTT'e 2a + F^RR' =A 2 N r , (3.20) 
-FT' 2 e 2a + F^R' 2 =A 2 (3.21) 



11 



and obtain 



—FTT'e 2a + F~ l RR' 
N r = =— =— 3.22 

e a (TR' — RT') 

N = - K J (3.23) 

V-FT' 2 e 2CT + F~ l R' 2 

As discussed by Kuchaf in Section IVA of [11], one can ensure that TR' — RT' and hence 
the Lapse function N are positive by an appropriate choice of x. y is then given from the 
definition (I3.12p as 

FT'e a 

(3.24) 



\/ —FT' 2 e 2a + F~ l R' 2 ' 
from which we obtain 

T'e CT = ^, (3.25) 
F 

where we used Eq. f)3.2ip . Using this to eliminate T'e a in Eq. (|3.2ip . we obtain 

F = -y 2 + !L. (3.26) 



as required by consistency with (13.101) 



In the above, we derived expressions for the generalized Schwarzschild time T in terms of 
the canonical ADM variables. As we will see in the following this determines the conjugate 
momentum to the Misner-Sharp mass function in a form that is appropriate for slicings that 
approach the Schwarzschild form at spatial infinity. Other asymptotic forms for the slicings 
are possible, including flat slice or generalized Painleve-Gullstrand (PG) coordinates: 

ds 2 {2) = -e 2a dT 2 G + (dR + Ge a dT PG ) 2 , (3.27) 

where a = a(T PG ,R) and G = G(T PG ,R). The geometrodynamical variables appropriate 
for such slicings were first derived in [37] . Since we have 

(DR) 2 = 1 - G 2 (3.28) 

for the above form of the metric, it follows that 

G = ±Vl — F. (3.29) 

By inspection of (I3.27p one can see that the positive sign yields an equation for ingoing 
null geodesies that is regular at any horizon F = 0, so this is the choice that is suitable 
for describing the spacetime near a future horizon (black hole). The opposite sign must be 
chosen for a past horizon (white hole) . We now go through exactly the same derivation as 
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before. Performing the coordinate transformations T PG = T PG (t,x) and R = R(t,x) in the 
metric (I3.27P and comparing to the ADM form (13.61) yields: 

A 2 = (B! + e a GT PG f - e 2a T PG , (3.30a) 
N 2 -A 2 N 2 = e 2a f 2 G -{R + e a Gf PG ) 2 , (3.30b) 
A 2 A" r = (R' + e a GT PG )(R + e a Gf PG ) -e 2u Tp G f PG . (3.30c) 

Solving (I3.30p for A" and N r , we find 

N = (R' + e°GT PG )(R + e a GT PG ) - e 2 °T PG f PG 
(R' + e°GT> G y-e°(T PG ) 2 



R'e a TpQ - Re a Tp_ 

\2 



N = = = = ^ PG =. (3.31b) 

\R' + e-GT PG f-e-{T PG ) 



To complete the derivation, we use ( 13 . 1 2 j) and the above expressions for A, A" and N r to 
calculate 

yA = {l-G 2 )e a T PG + GR', (3.32) 

which yields: 

e°T PG = y —± — ^ . (3.33) 

The second term on the right-hand side of the above guarantees that the PG time is well 
defined either for (with a +ve sign) future or (with a -ve sign) past horizons. 



4 Canonical formalism in general relativity 



In this section, we perform the canonical analysis for spherically symmetric spacetimes 
(k = 1) in general relativity without a cosmological constant, which is a generalization of 
the Kuchaf's analysis in four dimensions to arbitrary dimensions. We set 5m = 1 in this 
section for simplicity. The reduced action f !3.14j) then becomes quite simple: 



'M(GR) 



A 



J d 2 y 2{n - 2)R n -*y(N' r A + N r A') - 2N ((R n ~ 2 )" A~ l + (iT" 2 )'^- 1 )'^ 



+ (n - 2)(n - 3)(1 + F)NAR n ~ - 2(n - 2)R n - i yA 
Here A n -2 is the surface area of an (n — 2)-dimensional unit sphere, namely 

2vr (--l)/2 



(4.1) 



A 



n-2 



r((n-l)/2) 



V. 



(1) 



n-2ti 



(4.2) 
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where T(x) is the Gamma function. The purpose of this section is to show that the areal 
radius and the Misner-Sharp mass are well-defined canonical variables in the system, which 
will be generalized in Lovelock gravity in the following section. 



4.1 ADM variables 



We first derive the expressions of Pa and Pr. The corresponding Lagrangian density of the 
action ( 14. ip is 



C 



(n - 2) An-2 
2kI 



2R y(N' r A + N r A') 



( (R n - 2 )"i 



+ (n - 3)(1 + F)NAR n ~ - 2R n ~ A yA 



(4.3) 



from which we obtain P/v = Pn t 
_ _ (n-2)A n - 2r>n 



and 

R n ~ 3 N~ l (R — N r R') 



(4.4) 



R 



(n - 2)^-2 



2 ft? 



2R n ~ 3 N- 1 (NlA + N r A') - 2{n - 3)N- 1 AR n ~ 4 (R - N r R') - 2P n " 3 A^ 1 A 

(4.5) 

With Pa and Pr, the Hamiltonian density %{:— APa + RPr — C) is given by 



H ={N r AP A )' - N r AP' A + N r R'P R - 

(n - 2) An-2 



Nt R n - 3 1 -R"A- L + P'A~ Z A' ] 4 



n - 2)A n - 2 R n - 2 

n — 3 



Pa ( RPr ~ — APa 



"A-l i D'A-2A' 



AP r 



1 - A~ 2 P' 



(4.6) 



The first term in Eq. (14. 6 p is the total derivative and becomes a boundary term. Since 
P/v = P/v r = 0, the Hamilton equations for N and N r give constraint equations H = 
and H r = 0, where the super-momentum H r (:= &H/dN r ) and the super-Hamiltonian 
H(: = dTL/dN) are given by 



H r = - APi + R'P, 



R, 



H = - 



K~ 



[n - 2)A n - 2 R n - 2 



PARP R -— — AP A 



(n 2) A-2 (_#/ A -i + # A -2 A A + !L^ A P"- 4 f 1 - A- 2 P' 2 



(4.7) 



(4i 
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The action is finally written as 

4t(GR) = dt dx(AP A + RPr 



— NH — N r H r ). 



(4.9) 



It can be verified that with suitable boundary conditions the constraints H and H r are 
first class in the Dirac sense and generate spacetime diffeomorphisms that preserve the 
spherically symmetric form of the metric. 



4.2 The Schwarzschild-Tangherlini spacetime in various coordi- 
nate systems 

In this subsection, we review various coordinate systems in the Schwarzschild-Tangherlini 
spacetime. The Schwarzschild-Tangherlini vacuum solution in the best-known Schwarzschild 
coordinates is given by 



d* = -[l- ^-)de "V + r\^dz\ (4.10) 



where M := 2K 2 n M/[(n - 2)A n - 2 ] and M is the ADM mass. The corresponding ADM 
variables are 

^ = 1-^3- * = <>, ^=(l~y\ H = r. (4.11) 

In the next subsection, we will consider the boundary condition at spacelike infinity with 
this slicing. However, there is a variety of slicings in the Schwarzschild-Tangherlini space- 
time, as presented below. 

By introducing a new spacelike coordinate p as r = p[l + M/(4p n " 3 )] 2 /( n - 3 ), the met- 
ric (I4.10P is transformed into the isotropic coordinates: 



ds 2 = ~ 5 1 M J {Apn ~'I de + 1 1 + -*t- ) (dp' + p> lah dz a dz\ (4.12) 



4/(n-3) 

[l + M/(4p"" 3 )] 2 ^ ' ^ ' 4p"" :i ' 
The corresponding ADM variables are 



... 4/ (n-3) / ~ \ 2/(n-3) 

[1 + M/(4p"" 3 )] 2 I 4p«- 3 / H \ 4p— ! 



(4.13) 
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On the other hand, by introducing a new time coordinate r defined by 



dr := dt + \l 1- , (4.14) 

V r n - 3 (1 - M/r"- 3 ) 

the metric f)4.10p is transformed into the Painleve-Gullstrand coordinates: 



(h- = - ( 1 - -^3) dr 2 + 2 J ^drdr + dr 2 + r 2 lab dz a dz b . (4.15) 



The corresponding ADM variables are 



iV 2 = l, N r = \j^L, A 2 = l, R = r. (4.16) 



By the coordinate transformation 

" \ Ml/(n-Df^_Zij (417) 

from the Painleve-Gullstrand coordinates (14.15p . we obtain the Schwarzschild-Tangherlini 
metric in the Lemaitre coordinates: 



ds * = -dr 2 + M 2 /^- 1 ) 



2(n-3)/(n-l) / . \ 4/ (n-1) 



-(p - r) ) dp 2 + ( ::i 7r ^(p - r) ) labdz^z' 

(4.18) 



In this coordinate system, the central curvature singularity and the black-hole event horizon 
are represented by r = p and 

2MV(n-3) 

p-r = — , (4.19) 

n — 1 

respectively. In the Lemaitre coordinates, the radial coordinate p does not coincide with 
R in the asymptotically flat region. Defining a new radial coordinate as 

/ _ i \ 2 /(™-i) 
f:=MVM) ^ , (4.20) 



which coincides with i? in the asymptotic region, we transform the metric ( I4.18P into the 
following form: 

, 2 l2 / (n-l)M 1 / 2 r\- 2(n - 3)/(n - 1) , 2 „/ (n-^M^rX 47 ^ , a , „ 

^ = " ^ + (l " ^ 2f (n-l)/2 J ^ + r " 2f (n-l)/2 j ^ ^ " 

(4.21) 
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This metric is inhomogeneous and time-dependent in this coordinate system and the cor- 
responding ADM variables and their fall-off rates are 

N 2 =1, N r = 0, (4.22) 
/ (n-l)J^M -^/M) / (w -i)Mi/2 r y/(n-D 

\ 2r(™- 1 )/ 2 I , n-r\i 2f{ n-i)/2 ) ■ 



Lastly, we present the Schwarzschild-Tangherlini metric in the Kerr-Schild coordinates: 



ds 2 = - di 2 + dr 2 + r 2 lah dz a dz h + -^{di + dr) 2 (4.24) 
1 ' dP + ^-didr + ( 1 + — — j dr 2 + r 2 j ab dz a dz b . (4.25) 



q^TL 3 / iyTl 3 V iyTl 3 



The corresponding ADM variables are 

N 2_f 1 ^y 1 M f MV 1 A 2_ , ^ jR _ r / 4 26 ) 



4.3 Boundary condition and boundary terms 

To perform the geometrodynamics, the boundary condition plays a crucial role. In the 
present paper, we adopt the following boundary condition at spacelike infinity x — > 

JV + 0(aT ei ), (4.27) 

AT r ~A r °°(t)x- (n - 3) / 2 - e2 , (4.28) 

A~l + Ai(t)aT (n - 3) , (4.29) 

R-x + Rx (t)x~ (n - 4) - £4 , (4.30) 

where e\ is a positive number and e2 and €4 satisfy £2 > max[0, — (n — 5)/2] and €4 > 
max[0, — (n — 5)]. (The validity of this boundary condition is verified in Appendix[B]) The 
asymptotic behavior of Pa and Pr are given by 

P A ~ - (n ~^- 2 iV 00 1 ^x 1 " 8 * - A~x^ 3 )/ 2 ^ , (4.31) 

n - 2)A n -2 



P R ~ - ^ i^^A, 



2 



W)(-M + £ 2 )x("- 5 )/ 2 - e2 +Ax(t) 



(4.32) 



1 This boundary conditions are suited to asymptotically Schwarzschild slicings. The analogous boundary 
conditions for PG coordinates are given in |38j and discussed in Appendix 151 



17 



Under the boundary condition adopted, the Misner-Sharp mass converges to a finite value 
M ~ M°°(t), where M°°(t) is related to Ai(t) as 



Ai(t) 



«jjM°°(t) 
(n-2)A n -2 



(4.33) 



Now let us consider the boundary term for the action (14. 9p . The role of the boundary term 
is to subtract the diverging at the boundary in the variation of the above action. The action 
is completed by adding the boundary term, which gives a finite value in the variation. 

Since the variation of /m(gr) gives 



51 



M(GR) 



dt dxi d t (5AP A ) - SAP A + A5P A + d t (SRP R ) - SRP R + R5P R 



- SNH - N5H - 5N r H r - N r 5H r , . 



(4.34) 



we need to know the contributions from N r 5H r and N5H. Using the following results; 
N r 5H r = - N r 5AP A - (N r A5P A y + (N r A)'5P A + (N r 8RP R )' - 5R(N r P R )' + N r R'5P R , 

(4.35) 

N6H = ^irrelevant terms^ - ^ - ^" 4 "~ 2 j - ^(^"~ 3 A' 1 ^) / - (NR^A-^'SRj 

+ (N'R^A^SR)' - (N'R n - 3 A- 1 )'5R+ {NR n - 3 R'A- 2 5A)' - {NR n - 3 R'A- 2 )'5A 

(4.36) 

we can write (I4.34p in the following form: 



51 



M(GR) 



dt J dx ^dynamical terms J + J dx 



5AP A + 5RP R 



t=t 2 



i=ii 



dt 



-N r A5P A + N r P R 5R 



(n - 2)An- 



-NR^A-^iR') + N'R^A^SR + NR^R'A^SA 



l rm-3 A -1; 



?n-3 -Dl A -2; 



(4.37) 



Now the boundary condition comes into play. We assume 5 A = 5R — at t — t\,ti 
and then the second term in the above variation vanishes. Using the boundary condition 
(I4.27p - fl4.32p . we can show that only the contribution in the last integral comes from 
NR n - 3 R'A- 2 5A as 



NR n ~ 3 R'A- 2 5A ~ JVootfAi 



nlN^M^jt) 
(n-2)An- 2 



(4.38) 
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where we used Eq. (I4.33p . Finally we obtain the boundary term in a simple form: 



5L 



M(GR) 



dt I dx ^dynamical terms^ + J 



dt 



N^SM 00 ^) 



(4.39) 



4.4 Misner- Sharp mass as canonical variable 



In the ADM coordinates, the canonical variables are {A, Pa; R, Pr\. However, the physical 
meanings of the variable A is not so clear. In this subsection, we show that the two- 
dimensional equivalent action is written in a rather elegant manner by introducing the 
Misner-Sharp mass M as a canonical variable. We introduce a new set of canonical variables 
{M,P M ;S,P S } defined by 



S :=R, 

P s :=Pr-^(AP' a + PmM') 



M :-- 
Pm ' z 



in 



R 

2)A»-5 



2^n 



R 



n—3 1 



T'e a = 



yA 
F ' 



(4.40) 
(4.41) 

(4.42) 
(4.43) 



where we used Eq. (I3.25p . We are going to show below that, under the boundary condition 
(I4.27p ~ (I4.32p . the transformation from a set of variables {A, Pa;R, Pr} to another set 
{M, Pm', S, Ps} is a well-defined canonical transformation. 



Note (14.431) chooses the conjugate to M in terms of the Schwarzschild time T. As verified in 
Appendix [B] this leads to a finite Liouville form providing one chooses boundary conditions 
such the metric approaches the vacuum Schwarzschild solution sufficiently rapidly at spatial 
infinity. In order to use asymptotically PG slices, it is necessary to choose the conjugate 
to M in terms of the PG time Tpq. That is: 



P 



M 



yA y/T^F 



FR' 



(4.44) 



Since the extra term on the right is just a function of R and M, this corresponds to a 
straightforward canonical transformation (M, Pm, S, Ps) — > (M, Pm, S, Ps)- For simplicity, 
we henceforth stick to the Schwarzschild expressions. 



From the expression (I4.42p for the Misner-Sharp mass, we obtain 



PmM 



(n - 2) A n -2 yA Dn 



2«S 



R" 



F- (n-3)(l-F) 



(4.45) 
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which shows P M M = P A A+(- • • )R+S(- 
in terms of {A, Pa; R, Pr} as 



) + ( 



)'. The Misner-Sharp mass M is expressed 



M 



2 p2 



+ 



(n - 2)A n - 2 D „ 



R' 2 



(4.46) 



2(n - 2)A- 2 P n - 3 ' 2< 
From this expression, we can show that M' is a linear combination of the constraints: 

M' = A -1 (y# r - P'iPj , (4.47) 

where we used Eq. (14.41) to replace Pa by y. This implies that in the vacuum theory M is 
a constant on the constraint surface, as expected. 

Also using the expression (I4.46p . we can show that two sets of variables {A, Pa; R, Pr} and 
{M, Pm] S, Ps} satisfy the following Liouville form: 



P A 8A + P R SR = P M 5M + P S 5S + 5rj + (', 



where 



r] :=AP A + 



jn - 2) A n-2 

2kI 



R n ~ 3 R' In 



R' + yA 



(n - 2) An-2 un 



2kI 



P n " d ln 



R'-yA 
R' + yA 



R'-yA 



5R. 



(4.48) 

(4.49) 
(4.50) 



Under the boundary condition (I4.27p -( l4.32p . the total derivative term ( converges to zero 
at spacelike infinity. Hence, the transformation from a set {A, Pa; R, Pr} to {M, Pm] S, Ps} 
is indeed a canonical transformation, namely 

/oo poo 
dx(P A 5A + P R 5R) - / dx(P M 6M + P S 5S) = 5u[A, P A , ; P, Pr], (4.51) 
-oo J —oo 

/oo 
dx V [A,P A ,;R,P R ] (4.52) 
-oo 

is satisfied. It is shown that the integrands in the above equation, namely P\<5A + Pr5R, 
PmSM + PsSS, and r\ converge to zero faster than 0(x _1 ) at spacelike infinity under 
the boundary condition we adopt, and hence the above expression is well-defined. (See 
Appendix IB1 for the proof.) 

We now derive the Hamiltonian constraint and the diffeomorphism (momentum) constraint 
in terms of the variables {M, P M ; S, P$}- A straightforward calculation using the above 
equations verifies the following relation; 



£ - P M M 



NA 

~r7 



M' = - §PmM' + (t.d.) 



:^(V r + ArJ0M'+(t.d. 



(4.53) 
(4.54) 
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where (t.d.) is a total derivative term, we obtain the Hamiltonian density Hq in the equiv- 
alent two-dimensional theory as 

U G :=P M M + P S S-C 

=N M M' + N S P S , (4.55) 

where we have used Eq. (14.541) and defined new Lagrange multipliers N M and N s as 

R 



^(N+^(Ny + N r R^ 



N[^-^f)+N r P M , (4.56) 



y R' R' 

N s :=S = (Ny + N r R') . (4.57) 



Collecting terms in N and N r , we can express the total Hamiltonian as 



N 



\(P M y 




[\ Rf 


R'J 



N r (P M M' + P s S'). 



(4.58) 



The coefficients of N and N r are the super-Hamiltonian H and the super- momentum H r , 
respectively. These can in principle be expressed in terms of Kuchaf ' variables by doing the 
inverse canonical transformation. Note that one can replace H by the linear combination 
of constraints 

y ~ A 



Q\=H- —H r = -— M' 
R' R! 



(4.59) 



in agreement with (I4.47p . 

The Lagrangian density for the canonical coordinates (M, S, N M , N s ) is now written as 

C = P M M + P S S - N M M' - N S P S} (4.60) 

which corresponds to Eq. (122) of JIT] . The constraints, M' = and Ps = 0, are obtained by 
varying the Lagrange multipliers iV M and N s , respectively. On the constraint surface M = 
m(t) and Ps = hold. The reduced phase space is therefore two-dimensional consisting of 
Pm '■= dxPM{x,t) and m. With suitable boundary conditions (39], one can repeat the 
analysis of [UJ for spacelike slicings that intersect both left and right branches of the outer 
horizons of eternal black holes to obtain the reduced action: 



'(2) 



dt 



Pm'fh — (N + — N-)m 



(4.61) 



where N± := =p \im x ^± OQ N M . The reduced equations of motion in vacuum then imply that 
m = m = constant, and p m = —(N + — iV_). 
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Lastly, let us derive the boundary term in Eq. (I4.39P with the new canonical variables. 
Starting from 



'M(GR) 



dt / dx(P M M + P S S - N M' - N b P s ), 



(4.62) 



we obtain 



5L 



M(GR) 



J dt J dx\SP M M + d t {P M SM) - P M SM + 5P S S + d t (P s 5S) - P S SS 
- 5N M M' - (N M 5M)' + N M '5M - 5N S P S - N s 5P s ^j 



dt J dx ^dynamical termsj + J dx 



P M SM + P S SS 



t=t 2 



t=ti 



- / dt 



N M SM 



(4.63) 



Under the boundary condition fl4T27l) -f l432l . we obtain 5M ~ 5M°°(t) and N M ~ -N^t) 
at spacelike infinity and hence we obtain the same result (14.391) by setting 5M = and 
8S — at t — ii, ti- One important advantage of the new set of canonical variables is to 
greatly simplify the calculations. We will take advantage of this simplification in the next 
section. 



5 Canonical formalism in Lovelock gravity 



In this section, we show that all the results in the previous section can be generalized to full 
Lovelock gravity. In particular the transformation from the ADM variables {A, Pa; R, Pr} 
to {M, P M ; S, Ps} is a well-defined canonical transformation using definitions of Pm, S, 
and Ps that are the same as those in general relativity, Eqs. f l4.40p - fl4.43p . and M defined 
by Eq. (J2ZU). 
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5.1 ADM variables 



First we derive the ADM conjugate momenta Pa and Pr. The Lagrangian density from 
the action (I3.14p is 



(n- 2)A n -2 v-^ _ 
2^ ^ a[p) 

n p=0 



2pR n - 2p - 1 y(N r A)' 



2pN 
n — 2p 



( R n-2py A -l 



+ (n - 2p - 1)|(1 - Ff + 2pF^NAR n - 2 ~ 2p - 2pR n - 2p - 1 yA 
+ pR n ' 2p - 1 {l - (1 - Fy- 1 \{(AN r y + A-'NR')^ - Ay^ 



(5.1) 



Using the binomial expansion and integration by parts many times, we can rewrite the 
above Lagrangian density into the following form up to the total derivative. The derivation 
is presented in Appendix IA.2I 



f o\ A [«/2] 

[n - 2)A n - 2 ~ 

c = — — — 2^ a (p) 



2k1 



2pR n - 2p - 1 y(N r A)' 



2pN 
n — 2p 



( jR n-2 P y A -l 



V 



+ (n - 2p - 1) { (1 - F) p + 2pF }NAR n ~ 2 ~ 2p - 2pR n - 2p ~ 1 yA 



(n - 2) A-2 



[ra/2] j-p-2 



p!(-l) 



p— 1— U) 



n p=2 L u;=0 v ^ y 



AN r yR n - 2p - l F p - 2 - w (A~ 2 R 



2 D' 2 \/ 



p—2—w 



+ E 



2(p - 2 - w)!(-l) p - 2 -"'^ (A^ rJ R™- 2 P- 1 A 1 -^ J R ,2j )V(P- u '-^- 1 



j!(p -2-w-j)\ 



3=0 

(A- 1 NR'R n - 2p - 1 )'F p ~ 1 - 



2(p-w- j) - 1 



p-i 



-E 



p — 1 — w 
2p\ 



p-i 



E 



2 P \{-iy 



' w\(p — 1 — w)\ 



R n - 2p - l F w - l yA- 2 AR' 2 



w—l 

E 



jj2tt-l-j r g f ( i? n-2p-l A l-2j) jR /2j y 2( U ,-j)+l 



. 2( "^ +1 (2tw - 2j)!(-l) 9 R w - 'j+^R^p-^A 1 - 2 ^- 2 ^- 1 ^- 1 ^^)' 
3 ^ q\{2w-2j + l-q)\ 

2{w-j)-\ 



E 

9=0 



2{w-j + l) 

(2w - 2j - l)\(-l) q £2w-2j-q( pn-2p-l a -l-2j 7y-2(«;-j)+l p/2j+l+g yqy 

g!(2w - 2j - g)! V ' w 



(5.2) 



23 



From this Lagrangian density, we obtain 

.[n/2] 



Pa =- 



(n - 2)An- 2 



P-1 



X 



E 



pK-i) 



p=0 



A 2 



p=2 

(x _ 2j)y 2 ^ R' 2j - 2 



' w\(p — 1 — w)! 



iu— 1 



j!(w-l-j)! 2(w-j) + l A 2i 



-2 



(5.3) 



and 



/ A [n/2] 

p=0 



(N r A)' - A 
iV 



+ (ra-2p-lU (1-F) 



7 ^ 



(n - 2) An-2 



[n/2] 



-p-2 



2kI 



a ip) p (P ..,Z )!( / )P ...', t " ^"M 2(A- 1 iV J R'P n -^ 1 ) / F p - 2 -^ 



p=2 



4o=0 



— 1 — «;)! 



p—2—w 



+ E 



If — I ^p-2-iu-j 



i=o 



2(p-2- iy)!(-l) 
j!(p - 2 - w - j)! 



(N r ir- 2p - 1 A 1 - 2i B! 2i )'y 2 <r- v '- :i - 1 '> 



+ AA^ r i? ri - 2p - 1 F p - 3 ' w (A~ 2 i?' 2 ) / ^F - 2(p - 2 - w)y 2 ^ | 

p— 1 

2(p . - i)i(-ir r ^-2p-i A ~2 A RI 2 N - lFW - 2 f p _ 2/ _ 1} 2 \ 



U>— 1 



V" v L) ( d t (IF'~ 2p ~ 1 A 1 ~ 2j ) 



2(w-j)+l 

- E 

(7=0 
2(w-j)-l 



i a i-2n R' 2j y 2{w - j) n-2p-l Rn~2p-2 R > 2 i y ^-3)+i 



2(w -j) + l A 2 !' 1 



2j(2w-2j)\(-iy 



R 2(w-j)+l^ R n~2p-l^l-2j N ~2(w-j)-l R /2j-l+q N qy 



- E 



q\(2w - 2j + 1 - g)! 

(2w - 2j - ^)K~^) Q R 2w~2j-q-l ^n-2p-l ^-l-2j j^-2(w-j)+l j^/2j+l+q j^ q y \ 



- q\(2w-2j-l-q)\ 



(5.4) 



In general relativity, Pa = Pa [A, y(R)] and Pr = Pr[A, y(R)] can be algebraically solved to 
give a unique set of A = A [Pa, Pr] and R = R[P\, Pr]- In higher-order Lovelock gravity, by 
contrast, it is not possible to obtain a unique expression in general because of the fact that 
Pa = P\ [A, y(R)] and Pr = Pr[A, y(R)] are higher-order polynomials of y. As a result, it is 
difficult to obtain the explicit forms of the super-momentum H r and the super-Hamiltonian 
H, such that 

C = AP A + RP R -NH - N r H r . (5.5) 



24 



in terms of the ADM variables. However, it is not necessary to do so at this stage. Things 
are greatly simplified by using the generalized Misner-Sharp mass as a new canonical vari- 
able. As we will show, the super-momentum and the super-Hamiltonian with the new set 
of canonical coordinates are the same as those in general relativity and then the boundary 
terms at spatial infinity can be easily derived. 



5.2 Generalized Misner-Sharp mass as canonical variable 



We introduce a new set of canonical variables {M, Pm\ S, Ps} defined in the same way as 
in general relativity, namely by Eq. (12. lip , and Eqs. (14.40}) — ( 14.43[) . Then, we prove that 
{A, Pa; R, Pr} and {M, Pm', S, Ps} again satisfy the Liouville form (14.481) with the following 
total variation and the total derivative terms. The derivation is presented in Appendix IA.3I 



rj :-- 



(n - 2)An- 2 


[n/2] 

a (p) pR n - 1 - 

L p =i 




[n/2] 


P-i 


- E & (p) Rn 

p=2 


-l-2p 

' w(p — 

w=l v 


(n - 2)An- 2 


An/2] 

E a (p) pR n - 1 - 


2kI 



R' + yA 



pi 



w—l 

E 



R'-yA 

2(-i) 2 ^- 1 -i_R /2j y 2 ("'^')+ 1 



p-i 



E 



-p=i 

2p\(-l) w yR' 
w)\A 



R' + yA 



R'-yA 



' w\(p 



w—l 



F w ~ l + 



j=0 



[n/2] 

+ e ^ P) R n - 1 - 2p 

p=2 

2j(w - l)\(-iy-^-jR> 2 i- 2 y ^-i) 
j\(w - 1 - j)\[2(w - j) + 1]A 2 ^'- 2 



(5.6) 



SR. (5.7) 



In Appendix \AA\ it is proven that Eq. (14. 54[) still holds in full Lovelock gravity. This 
immediately implies that the Hamiltonian density in the equivalent two-dimensional theory 
takes the same form as that in general relativity (I4.55p . where the definitions of the new 
Lagrange multipliers N M and N s are the same as those in general relativity (I4.56P and 
(I4.5?p . Finally, the Lagrangian density for the canonical coordinates {M, P M ; S, P$} can 
be again written as 

C = P M M + P S S - N M M' - N S P S (5.8) 
and the super-Hamiltonian and super-momentum constraints are again as in (14.581) . 

In comparison to the rather complicated starting point in Eq. (13.41) . this equivalent La- 
grangian density is extremely simple and the physical meaning of the canonical variables 
are very clear. Remarkably, the coupling constants do not appear explicitly in (15. 8p . 
They are in fact hidden in the definition of the mass function. This makes it possible to 
treat any class of Lovelock gravity in exactly the same way. 
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5.3 Fall-off rate at infinity and boundary terms 



In order to prove that the transformation from {A, Pa, ; R, Pr} to {M, Pm\ S, Ps} is canon- 
ical and well-defined, we have to discuss the asymptotic behaviour of the variables. We 
adopt the same boundary conditions fl4.27p - fl4.30p as in general relativity. With these con- 
ditions, one can verify that the generalized Misner-Sharp mass (12.111) behaves near spacelike 
infinity as 



M 



(n - 2)A.-2«(i)Ai(t) 



ft", 



(5.9) 



This is the same as in general relativity and hence we set Ai as in Eq. (I4.33P (where a 
in order that M ~ M°°(t) at infinity. 



(i) 



It can then be shown that the leading terms of Pa, Pr, C, Vi Ps, Pm, N M , and N s are the 
same as those in the general relativistic case under the boundary condition (I4.27p - p.30p . 
As a consequence, the proof carries over from general relativity and all the terms in the 
Liouville form (I4.48P are well behaved at spacelike infinity. This is sufficient to prove the 
transformation from {A, Pa, ; R, Pr} to {M, Pm', S, Ps} is indeed a well-defined canonical 
transformation, and that the Hamiltonian dx(N M M' + N s Ps) is also finite. 



We now have the following two-dimensional action with a new set of canonical variables; 



'M(L) 



dt 



dx{P M M + P S S - N M M' - N S P S ) 



(5.10) 



with the same asymptotic behavior as in general relativity. The boundary term for the 
above action that makes the variational principle well defined is then also the same as in 
general relativity: 
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M(L) 



dt 



dx ^dynamical terms^ + J 



dt 



N^SM^it) 



(5.11) 



Given the above, we can now write down the super-Hamiltonian and super-momentum 
constraints for full Lovelock gravity. In terms of the geometrodynamical variables they are 
the same expressions as in general relativity: 

H = {jf-^)M' + yP s , (5.12) 

H r = P M M' + P S S'. (5.13) 

The expressions in terms of ADM variables are considerably more complicated and can in 
principle be obtained once again by substution from Eqs. (I4.40p - (l4.43p . with M defined by 
Eq. d2HU>. 
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6 Adding matter fields 



In this section, we introduce matter fields in the argument with the ADM variables dis- 
cussed in the previous sections. Here we write super-momentum and super-Hamiltonian 
for gravity as and in order to distinguish from the total super-momentum and 

super-Hamiltonian including matter contributions. The following argument is valid in full 
Lovelock gravity. 

It can be shown from Eqs. (|512j) . fOTj) . (jOUjh ffl~43l and ffl~54|) that the gravitational 
Hamiltonian Hq is given by: 

H G = J dx(NH<® + N r H^), (6.1) 

where 

HP =P s S' + P M M' = PrB! - P' A A, (6.2) 

tf(G) = _A M ' + |-#( G ). (6.3) 
HI K 

We have used (I4.47|) to derive (16. 3p . Since Eq. (15. 3p shows that y is not a function of N 
or N r , we can see that the Hamiltonian density is the sum of Lagrange multiplier times 
constraints. 



6.1 Massless scalar field 



First we consider a massless scalar field if) as a matter field, of which action is /matter = Itp 
in the action 02. ip : 

h = ~\\ d n x^(Vif>) 2 . (6.4) 

The equivalent two-dimensional action in the symmetric spacetime under consideration is 
given by 

h = J d 2 y^g^R^(D^ (6.5) 

This gives the momentum conjugate P^ to if) as 

P* = An - 2 ^ Rn ~ 2 (if) - N r if>') , (6.6) 
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with which we can write the matter action as 

1 



dxdtipP^ — / dxdtN 



2A 



p2 



An- 2 R n - 2 



+ A 



n-2 



R n ~ 2 ^ 2 



N 



(6.7) 



Equation (I6.5P does not contain any derivatives of the metric or P, which means that 
adding the scalar action to the gravitational action (14. 55ft does not change Pa or Pr. This 
allows us to write the total Hamiltonian as the sum of the gravitational and matter parts. 
Using Eqs. (16. ip . (16. 2p . (16. 3p . and (16 .7p . we obtain the total Hamiltonian P to tai as 



total 



dxN 



A M'+ V H + Nr H + 1 ( P l 



A- 2 P n ~y 2 ) + P^'§ 



dx 



N(H™ + P (M) ) + N r (HP + HM) 



(6.8) 



where y is a function of the phase space variables, A, Pa and R via Eq. (I5.3p . The super- 
Hamiltonian p( M ) and super-momentum for ip are given by 



2A \An- 2 R" 2 



p2 



+ A- 2 P"" V 2 



(6.9) 
(6.10) 



The Poisson bracket of Hamiltonian constraint, H = with the total momentum 

constraint, H r = HP + is given by 

{H, H r } = {H^ G \HP} + {P(°), Ff>} + {ffM H (G) } + {H M t H (M)y (6 n) 
Because our theory is diffeomorphism invariant, {H,H r } must be weekly equal to zero. 

Because there are two first class constraints, there are two gauge choices to pick. We choose 
our first gauge as 

X :=R-x^0. (6.12) 

This forces the spatial coordinate to be the areal radius which means that R is no longer a 
phase space variable, it is now a coordinate. In order to insist that x is satisfied at every 
time slice, we must insist that x — {x> H} ~ 0, which shows N r /N + y/R' « 0. We use 
this relation to write one Lagrange multiplier in terms of the other. This leaves us with one 
Lagrange multiplier which reflects the fact that there is only one gauge fix left to choose. 



We can now plug the gauge choice (16.121) and its consistency condition into the Hamiltonian 
as long as we use Dirac brackets to evaluate the equations of motion in the end. Note that 
the remaining phase space variables, A, P\, ip and P^, all commute with x an d so the 
Poisson bracket is the same as the Dirac bracket. Plugging x — X — Eq. (16.81) gives 



P, 



total 



dRN 



-AM' + 



p2 



2 A \A n -2R n - 2 



+A- 2 p n -y 2 



(6.13) 
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In the last term we replaced N r /N by — y as required. Since the mass equation (12. lip is 
written as 

M = (n = 2) 2 A - 2 £ a w JT" 1 -* (1 - A" 2 + y 2 ) P , (6.14) 

K ™ p=0 

we can write ?/(= N r /N) in terms of the mass function. For this reason we leave the factor 
of N r /N in the Hamiltonian with the understanding that it is the solution to Eq. ( 16.141) . 

For our second gauge choice we choose 

£:=A-1«0. (6.15) 

By the same reasoning used for the first gauge choice we can set £ strongly to zero (namely 
since A commutes with ip and P^) which gives the Hamiltonian 



tftotal = / dRN 

and the mass function 



M = ^^ n — Z^w* [W ■ (6 ' 17) 

To see the significance of this gauge choice, notice from Eq. (12. lip that g 11 — > 1— 2/t^M/[(n— 
2)*4.( n _2)tt(i)-R n-3 ] in the general relativistic case when we strongly set £ and x to zero. This 
gives the metric in the non-static version of Painleve-Gullstrand coordinates: 



(6.16) 



ds 2 = -N 2 (l- —) dt 2 + 2N* / 7 —-dtdR + dR 2 . 

(6.18) 

To ensure that the second gauge condition is conserved in time we must insist that d(A — 
l)/dt = {A — 1, H} = — > 5H/5P\ = 0. Although we have chosen to write N r /N in terms 
of the mass function, it can also be written in terms of Pa- All of the Pa dependence in 
the Hamiltonian is in the terms of N r /N. Therefore we can write 



<5Ptotal 5 



SPa SP a 



J dR \N'M + NP^'^A 



JM d(N/N) djNr/N) 
d(N r /N) dP A ^ dP A 



from which the consistency condition is given as 
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where it is understood that we use (16. 17)) to find dM/d(N r /N) and write N r /N in terms of 
the mass function M. Notice that the actual relation between N r /N and Pa is not needed. 



Using Hamilton's equations and Eq. ( 16. 16f) . we find 



ip =N 



A n -2R 



n-2 



Nr 

N 



N ( An^R n - V + 



(6.21) 
(6.22) 



These equations, along with the consistency conditions (16.171) and (16.20)) and the Hamilto- 
nian constraint 



— M' 



1 



p2 



2 



A n ^R n - 2 ^ + P^ 



N 



0. 



(6.23) 



determine the evolution of a collapsing scalar field. 



6.2 Charged scalar field 



In this subsection, we consider a U(l) gauge field coupled to a charged complex massless 
scalar field ip = (ipi + iip2)/V^i where ipi and ip2 are real functions. We write the action 
for this matter as / ma tter = Iem'- 



'EM 



d n xyf^g 



- (d» + ieA") r {d„ - ieAJ </> - -F^F^ 



(6.24) 



where e is the charge and the Faraday tensor F^ u is defined in terms of the gauge field as 
Fpv = d fl A u — d u A IM . Under the symmetry assumption in the present paper both for gravity 
and matter, the equivalent two-dimensional action is given by 



Iem = A n -2 J dtdx^-g {2 )R n 2 



(d B + ieA B ) r (d B - ieA B ) $ - -F AB F AB 



. (6.25) 



Adopting the ADM coordinates, we obtain 



Iem =—^— / dtdx — — — 



+ ^l) - 2N r (M[ + ^2) + (N 2 - iV 2 A- 2 )(^ 2 + 



- 2e|(A - AUiXV^i - M 2 ) - [ N r (A - N r A x ) + N 2 A~ 2 A 1 ) (#i - *p'M 



+ e 2 [ (Ao - N r A,) 2 - N 2 A~ 2 Af ) + ^) + A^A, - A' ) 2 



(6.26) 
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where A n dx^ = A Q (t,x)dt + A\{t, x)dx. From the above action we find the conjugate 
momenta: 



P 



01 



ip2 



AO 



Pax 



A n - 2 R n - 2 A 
N 



N 



i> x - N r ^[ + e(A - N r A{)^ 2 
i) 2 - N r ip2 -e{Ao- N r Ai)i/)i 



--0. 



•A n - 2 R n ~ (A x - A' Q 
NA 



(6.27) 

(6.28) 
(6.29) 

(6.30) 



which give the Hamiltonian for the present matter field 

N 



H 



EM 



dx 



(Pii + Ph) + e(^o - NrA^P^ - P^ 2 ) 



+ N r {P^[ + P^' 2 ) + 



NR 



n-2 



2A- 2 A 



(eAi^i-^2) +(eAi^ 2 + ^i 



l\2 



+ 



NA 



2A n - 2 R- 2P * 1 + PmA '° 



(6.31) 



Since the action (I6.26P contains no derivatives of the metric or R, the addition of Jem to 
the gravitational action does not alter the Hamiltonian analysis and allows us to write the 
total Hamiltonian as 



H 



total 



dx 



N(H™ + # (EM) ) + N r (H^ + Iff M )) + AqH 



(EM) 
AO 



(6.32) 



where H^ G \ H^ G) and Hf M) are given by Eqs. (J6l| and flOJ and H^ and ff r (EM) are 
given by 



h (em) 



p2 , p2 

'2A n - 2 AR 



n-2 



A n - 2 R n 
2A 



/ \2 



API 



p(EM) = _ ^(p ^ _ p ^ + (p ^ + p ^ 



(EM) 
AO 



=e(P^Vi - ^1^2) - ^ 



2A- 2 i? n - 2 ' 

(6.33) 
(6.34) 
(6.35) 



where we used integration by parts and asymptotic condition PaiA — > at infinity to 
derive Eq. f l6.35[) . The consistency condition on the constraint (16.291) is {Pao, H total } = 
0, which gives e(P^ 2 ipi — P^iip 2 ) — P' A1 = 0. This condition is already added into the 
Hamiltonian with Aq as its Lagrange multiplier. Since Pao is weekly equal to zero we can 
use the equation of motion for Pao to show that H A E Q M ^ is weekly equal to zero and is, 
therefore, a constraint in the same way as the constraints multiplying N and N r . 
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This Hamiltonian is composed of three first class constraints which means that there are 
three gauge choices to make. Our first two gauge choices will be the same as in section I6TT1 
Using similar reasoning we can write the Hamiltonian as 



Ptotal = / dR 



p2 _|_ p2 a t}n-2 



p2 



2A n - 2 R n ~ 2 



^ f-eA^P^i - ^1^2) + (P^i + ^2^2) J } 



(6.36) 



Just as in section 16.11 the consistency condition on the first gauge fix requires us to write 
N r /N as a function of M using Eq. (16.171) . The consistency condition on the second gauge 
choice, analogous to Eq. (I6.20p . is given by 

N ' d(N M /N) + N (- e MP^i ~ PtpM + (P^'i + Pwlti) = 0- (6-37) 
For our third gauge we choose 

e := At « 0, (6.38) 

which is the coulomb gauge with the constant, A\ chosen to be zero. This condition, along 
with the electromagnetic constraint, 

e(P^i - P^ 2 ) - P' A1 « 0, (6.39) 

removes A\ and its conjugate momentum Pa\ from the set of phase space variables. We 
can therefore set e strongly to zero in the Hamiltonian as we did for the first two gauge 
choices. This gives the following Hamiltonian: 



#total — I dR 



N 1 _ M , + p'i + g, + A- 2 ir- 2 ( ^ 2 + ^ 



+ 



p2 



2A n - 2 R 



■n-2 



+ ^(P^i + P^ 2 )}+^o(e(P^i-P^2)-PAi) , (6-40) 



Eq. (I6.38P is given by 



where it is understood that Pai is the solution of Eq. (16. 39 p . The consistency condition on 



NP 



eiV /" 

4. « - An 2jRn _ 2 y dfl(^ - ^2), (6.41; 



which puts a condition on the final Lagrange multiplier and must be satisfied at every time 
slice. This is the last consistency condition on e. 
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With the fully gauge fixed Hamiltonian (I6.40p we may write down Hamilton's equations of 
motion in terms of the remaining phase space variables, ipi, P^\, ip 2 and P^ 2 . The equations 
of motion are given by 



iPt —N 

ik =n 

Pip I 



A n . 2 R n ~ 2 N 



+ — f^i ) - eA ^2, 



P 



'tj'2 



A n - 2 R n ~ 2 N 



N 

+ ^2 ) + e^o^i, 



N r 



eA P r 



■02 5 



1p 2 



N [ A n . 2 R n - 2 ^ 2 + 



N r 
~N' 



+ eA Q P t 



(6.42) 
(6.43) 
(6.44) 
(6.45) 



It must be remembered that at every time slice the equations of motion must be supple- 
mented by the consistency conditions (I6.17p . (I6.37p . and (I6.4ip . as well as the Hamiltonian 
constraint: 



-M' + 



]^k + + + = o, (6.46) 



where it is understood that Pax is the solution of Eq. ( I6.39P . 



7 Conclusions 



In this paper, we have performed the Hamiltonian analysis for spherically symmetric space- 
times in general Lovelock gravity in arbitrary dimensions. We have shown that, as in general 
relativity, the areal radius and the generalized Misner-Sharp quasi-local mass M are natural 
canonical valuables that yield the remarkably simple, geometrical action (I5.8P for the generic 
theory. Using these variables also enabled us to rigorously derive the super-Hamiltonian 
and super-momentum constraints (I5.12p and (I5.13P for the most general theory, a task that 
would have be daunting at best, if not impossible, in terms of ADM variables. Most im- 
portantly, our results are useful: the geometrodynamic variables allow the physical phase 
space of the vacuum theory to be explicitly parametrized in terms of the ADM mass and its 
conjugate momentum, as done for general relativity by Kuchaf [II] . This in turn provides 
a rigorous starting point for the quantization of Lovelock black holes using the techniques 
of [12]. Finally, the simple form of the Hamiltonian allows us to gauge fix and derive the 
Hamiltonian equations of motion for the collapse of self gravitating matter in flat slice 
coordinates. We are now in a position to study the dynamics of black hole formation in 
generic Lovelock gravity. These investigations are currently in progress. 

Acknowledgments HM thanks Jorge Zanelli for discussions. This research was supported 
in part by the Natural Sciences and Engineering Research Council of Canada and by the 



33 



JSPS Grant-in-Aid for Scientific Research (A) (22244030). This work has been partially 
funded by the Fondecyt grants 1100328, 1100755 and by the Conicyt grant "Southern 
Theoretical Physics Laboratory" ACT-91. The Centro de Estudios Cientificos (CECs) 
is funded by the Chilean Government through the Centers of Excellence Base Financing 
Program of Conicyt. GK and TT are very grateful to CECs, Chile for its kind hospitality 
during the initial stages of this work. HM would like to thank University of Winnipeg, 
for its hospitality while part of this work was carried out. GK also thanks Jorma Louko 
for introducing him to the beauty of Kuchar' geometrodynamics. TT and GK also thank 
Danielle Leonard and Robert Mann for the collaboration in [26] which was instrumental in 
leading to the present work. 



A Derivations 

In this appendix, we present the details of several lengthy derivations of results in the main 
text. The following relations will be useful for much of the following: 




(A.l) 



(A.2) 



Equation flA.2j) is shown as 




■yA 



5((R' 2 -y 2 A 2 )/A 2 ) 
(R' 2 -y 2 A 2 )/A 2 




2y5A + yA 



R' 2 - y 2 A 2 ' 
2R'5R' 



2R'5R' 



2<%A) - 2y5A - 2R 



5(yA) 



yA R' 2 -y 2 A 2 



(yA) 2 




(A.3) 



Equations (IA.1I) and ( 1A.2I) will also be used by replacing 8 by d t or d : 
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A.l Lagrangian density (13.51) 



We now derive the Lagrangian density (13. 5ft from Eq. ( 13 Ah . Using the binomial expansion 
for the last two terms in (\3Ah yields 



2p(p-l) 

__ R -2( P -1) 



(D 2 R) 2 - (D A D B R) (D B D A R) (k- (DR) 2 \ p ' 2 (?) ( k - {DRf\ p ^ 



R 2 



R 2 



R 2 



i=0 

(D A R) (D>D A R - D A D 2 R) g 



(i + l)!(p-2-i)l 



i=0 



(i + l)!(p-2-i)! 



(A.4) 



where we used the following two-dimensional identity: 

(2) 

{DRfTl = 2{D A R) (D 2 D A R - D A D 2 R) . 
This identity can be derived from Eq. (2.10) of [26] . 

Eq. (13 .4p now reduces to 
n- 2) 



(A.5) 



(p) 



(n-2p)! 



(n — 2p){n — 2p — 1) 



k-(DR) 2 \ p n/ n ,D 2 R fk-{DRf^ p - 1 
1 —2p(n — 2p)- ' 



R 2 



R 



R 2 



(2) 



P-2 



+ pk v ~ x R 2 - 2v ll + ^ 



i=0 



2(-l) i A*- 2 -'p!(.D.R) 
(i + l)!(p-2-i)! 



2/ 



(D 2 R) 2 - (D A D B R)(D B D A R) 



R 2 p~ 



D A R(D 2 D A R - D A D 2 R) 1 



(A.6) 



Using integration by parts, we can rewrite the term in curly brackets in (1A.6j) as 



P-2 

E 

i=0 



- D A R(D 2 D A R - D A D 2 R) 



E 

i=0 



2(-1)^p- 2 -V! 
(i + l)!(p-2-i)!' 



(Di?) 2 ^^™- 2 ^ -D A ({DR) 2 ) - (D 2 R)(D A R) U(t.d.) 



(A.7) 



where we used the following identity: 



- (^b*)(^x*) + - ^Wgg s 0. (A, 



2(Di2) 1 



(£>J2)* 
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Using the above result together with the integration by parts and the following identity; 



|4 2(-l) i kP- 2 ~ i p\(DR) 2i 



i=0 



(i+l)!(p-2-i) 



£ 

w=l 



2(-l) w k p ~ l - w p\(DR) 2w (DR) 
w\(p — 1 — w)\ 



w=0 ' ^ >' 



2pk p ~ 1 - 2p(k - (DR) 



2\p-l 



(DR) 2 

we can rewrite expression (1A.6|) in the form (13. 5 p up to a total derivative. 



(A.9) 



A. 2 Lagrangian density (15.21) 



In this appendix, we show how to derive the Lagrangian density (15.21) from the action 
(I3.14p . While we consider only the spherically symmetric case (k = 1) in the main text, 
here we derive the equations for general k. 



For this purpose, we separate the action ( I3.14p into two portions: 
J M =/i + /2 + (t.d.), 



(A.10) 



(n 2W {k) [n/2] i 



2Ki 



p=0 



yoi(p) 



2pkP- 1 R n - 2 P- 1 y(N r A)' - 2pkP lN ( (R^YA' 1 ) 

n — 2p \ J 



+ pR n - 2p - 1 ^k p - 1 - (k - FY" 1 \(AN r y + A^NR')^ 
+ (n - 2p - 1) j(A; - F f + 2pk p - 1 F^NAR 



yn-2-2p 



(A.n; 



(n - 2)V {k) [n/2] r 

h ■■= — — — 2^ d yp a (p) R 



2Ki 



p=0 



-2F _1 yA - {jfeP- 1 - (k - Ff-^Ay- 



(A.12) 



In order to perform the variation of Im, we have to deal with the terms containing F' and 
F. Using the binomial expansion and integration by parts, we can rewrite the second line 
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of Eq. (TATTD as 



_ _ o\v W [n/2] ( 1 

K n p=Q I J 



2«* ^ ^ ^ «,!(p -!-«,)! ^ F 

MjxSKL ^ (p-l)!(-ir^ 
24 l^ a WPZ^ ^(p-i-u,)! 

" p=2 «)=0 v ^ ' 

P ^ W 2(p - 2 - w)\(-l) p - 2 - w -i (N r R n - 2p - 1 A 1 -^R' 23 )'y 2( ~P~ w -^- 1 
3 ] -{P ~ 2 - w ~ j)- 2(p-w-j)-l 

+ \N r yR n - 2p - l F p - 2 - w {A- 2 R' 2 )' - (hT 1 N R' R n ~ 2p ~ 1 )' \. (A.13) 

p — 1 — w J 

This expression does not contain y' and can be used to obtain Pr. An important observation 
is that the only part of the action Im in (IA.10I) that contributes to Pa is After a tedious 
but straightforward calculation, using binomial expansion and integration by parts yet 
again, we can rewrite the last term in I2 as 

(» - 2)V« '"' 21 



" p=0 



r ( 1 f 

J d'ypa^R^- 1 ^- 1 — (k — F) p ~ l jAy— 



in-2) V ^ fj2 _ „ p ^(p-m-irk p 



E / ^ E (p -j£lZ% 

p=2 ^ tu=l 7 ' 

(n-2)KS^ ] /• ^ (p - l)!(-l)»A*-i-» ^ ( W - 1)!(-1)"-W 

+ .2 2^ / a f^to w! ( p _!_„,)! 2^ j\(w-l- j)\ 



K ' n p=2 



x 



1 fUF?"-^- 1 A 1_2 ^ P /2 - ? „ 2 ( u '-j)+ 1 

2(w-j) + l 1 



^ Q q\(2w-2j + l-q)\ w-j + l K ' rJ 

_ ( (2w - 2j - l)!(-l) g R 2w -^-i 2p _ lA _ x _ 2j . 2(M _ 3 - )+1 Rl2 j+i+ q N gy 

^ g!(2 W -2j-l-g)!2^-2j-g l rJ 

+ $(•••) + 3*0 ••)■ (A.14) 

Using Eqs. (1A. 13h and (1A.14I) . we obtain the Lagrangian density (15.21) . 
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A. 3 Liouville form (14.481) in Lovelock gravity 



In this appendix, we verify the Liouville form (14.481) in Lovelock gravity. Note that the 
explicit form of Pr is not used in the derivation. 



Using Eqs. (EOT]) and (EH}, we write P M bM as 

[n/2] 



P M 8M 



p=0 

[n/2] 



[n 



1 - 2p)(l - F) P 5R - pR(l - F) P ^5F 



~ (n z ] / n ~ 2 E «W(« - 1 - 2p)^iT- 2 - 2p (l - F) P 5R 

n p=0 

I (n ~ 2).A w -2 - 3 2/A 

+ ^2 



+ 



(n - 2) A n -2 
2< 



p=i 

[n/2; 



P-1 



£ a^pir- 1 -* ^ i/ ;„ l) > r -//.\ / 



p=2 



-J w\(p — 1 — w)! 



(A.15) 



Using Eq. (|A.2[) for the second term and Eq. flA.ll) for the last term together with the 
binomial expansion, we obtain 



P M 5M 



[n - 2) An-2 
2 k 2 



[n/2] 

E 

p=0 

[n/2] 



«(p)l« 



1 - 2p)^-R n - 2 ~ 2p (l - F) P 5R 



+ 



+ 



(n-2)^_ 2 fe3 2 L D n-i-2 5 ? (P-I)K-Ir 



R' + yA 



R'-yA 



2kI 



p=2 L ui=l v ' 



w — 1 

E 

j=0 



- 1 



-y 



2(w-l-j) 



A 2 



-2^y)-2y(y 2 + F)5Aj 
(A.16) 



An important fact is that P^SM has the form of PmSM = P\5A + (■ ■ ■ )5-R + 5?y + £', where 
Sn and directly appear in the Liouville form (14 .48 p . because S is defined by S := R. 
Therefore all the terms with SA in Eq. (A.16) are contained in the expression for Pa- 
using the binomial expansion and integration by parts, we can calculate the other terms 
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in Eq. (1A. 16|) as follows 



[n/2] 



+ 



(n - 2) A-2 
2« 2 

(n - 2) A.-2 



p=0 

[n/2] 



I"/ -"J f 

E a^piT -1 -*] 2c%A) - #5 In 
p=i ^ 



R' + yA 



R'-yA 



[n/2] 



p-1 



(p-i)K-ir 



2 ,.2 p 1 ""E „• .. | 

" p=2 L ui=l v ' 



w—1 

E 

i=o 



(w-l)!(-l) 

- 1 - j)! 



-3/ 



2(to-l-j) 



R_ 

A 2 " 



/2 \ J 



yA 



5(R 



i2\ 



A 2 



- 2y5y 



[n/2] 



- ( '" 2 2 ^"' 2 E 5 <*><" - 1 - 2 p»f RMp{1 ~ F) ' SR 

n p=0 
[n/2] 

E a (p) 6 



+ 



(n - 2) An-2 
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- 2) An-2 



2kI 



p=i 

[n/2] 

E 

p=i 



n-l-2p 



2yA - i?' In 



R' + yA 



R'-yA 



2(n - 1 - 2p)R n - 2 - 2p yASR 



- in 



1 - 2p)R n - 2 - 2p R'5R + R* 1 - 1 - 2 ^^')) In 



R' + yA 



+ 



(n- 2)^-2^' „ D „-i- 2p v^ b-l)K-ir ^-iV* 



/ ■ 7 



[n/2] p-1 
p=2 ui=l v ' 



R'-yA 
F 



A 



5(R') 



n p=2 



' w!(p — 1 — io)! 



- 



i=o 



jl(w-l-j')! \A 2 J 2{w-j) + l_ 



, (n-2)A-2 ^L py 2 ^)+ 1 

f 5 — E a &» 2(ti; _ j) + i* 



K 2 , 



X 



tu— 1 

E 

J=0 



p=2 



p-i 



R n-l-2p A 



(p -i)i(-ir 

^ w!(p — 1 — ty)! 



VA 2 



(A.17) 



Because there will not appear any more total variation terms, we can read off the total 
variation term rj as Eq. (15.61) . 
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In order to derive £, we write down the quantity IT := P^SM — P^dA — Srj: 



n 



+ 



+ 



+ 



p=0 

[n/2] 



n - 2)Ai-2 



2k? 



n - 2)A n -2 



2< 
n - 2)^-2 



v 

[n/2] 

p=l 

[n/2] 

P =i 

[n/2] 



R' + yA 



R'-yA 



5R 



pR n-i-2p ln 



pR n-l-2p ln 




p-1 



n - 2)^4 n _ 2 



£a w 2 P g--*^ '?,- 1) ! ( - 1) '. F-gg« f 



p=2 
[n/2] 



P-1 



^ w\(j> — 1 — «;)! 



2kI 



(p -!)!(-!)' 



p=2 

[n/2] 



d w\(p — 1 — w)\ 



A 
A 



5R 



p-i 



(p -!)!(-!)' 



p=2 w=l K ^ ' 



X 



w— 1 

E 

i=o 



j!(w-l-j)! V A2 /2(u;-j) + l 



/v- 



[n/2] p-1 

y&MpR^-^Y , 



(«-2R-2^, .on-l-2p^ (P-I)K-I) 



X 



id— 1 

E 

i=o 



(to - 2jR' 2j ~ 1 A 1 - 2 i 
j\(w-l-j)\ 2(w-j) + l' 



y 



S(R'). 



(A.18) 



The last term generates a total derivative term by integration by parts. Now we see all the 
total derivative terms and can read off the total variation term ( to be Eq. ( 15. 7(1 . 

In order to prove the Liouville form ( I4.48p . we write down the quantity ESR := PmSM — 
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P A S A -Stj-C' as 

E5R = - (n ~ 2) 2 X ~ 2 £ a w (n - 1 - 2 P )^P- 2 - 2 »(l - F) P 5R 



*/2] 



™ p=0 
n/2] 



(w - 2) A-2 
2« 2 

(n - 2) A-2 



P =i 

[n/2] 



2(n - 1 - 2p)yA + P^ln 



R' + yA 



R'-yA 



p-i 



(p-m-iy 



10— 1 

E 

j=0 



p=2 u>=l v ^ 7 

I — r I — : SR 



j\(w-l-j)\ \A 2 J2(w-j) + l 



MAid 'y. Lfi.n-i-2^ (p-iJK-i) 1 
12 2^ a (p) PjK 2^ 



p=2 
ru— 1 



^ «;!(p — 1 — w)\ 



x 4 F 



+ E 



A 1 j\( w -i-j)\ 2 (w-j) + l 



-y 



2(w-i)+l 



5P 



(5/?. (A.19) 



We now show that 

S=i(A^ + P M M / ), (A.20) 
which is sufficient to verify the Liouville form (14. 48 p . From the expression (15. 3p . we obtain 



API 



- {H 2 \ An - 2 E a ip) pA{(n - 1 - 2p)R^yR' + BT^A 

K n p=1 I ) 

I _ 0\ A [n/2] 

+ ^ 2 J 2 E fi WP(n - 1 - 2p)P"- 2 - 2 *P'A 



x 



p=2 
p-1 



— 1 — w)\ 



^w!(p-l-w)!^ j\(w - 1 - j)\ 2(w - j) + 1 vy ; 



[n/2] 



2/€ « 2 



p-1 



' u>!(p 



"~ L (p-l)!(-l)« ^(u,-!)^-!)"-!-^ 2(1-2.7) 



^ w\(p -l-w)\^ 

w=l v ^ ' j=0 



^w\{p - I - w)\ j\(w-l-j)\ 2(w-j) + l 



y 



2(w-j)+l/ y 2 + py 



(A.21) 
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Using 



R'< In 



R' + yA 



2{yA)'-2yA'-yAy 



R'-yA 

for the logarithmic term, we finally obtain 

( -0\A [n/2] A 

E5R = - {H £ - 1 - 2p)^- 2 - 2p (l - FY6R 



(A.22) 



+ 



2< 
(n - 2) A»-2 

(n - 2)A.-2 



p=0 
[n/2] 

p=l 
[n/2] 



2p(n - 1 - 2p)R n - 2 - 2p yA + pR n ~ l - 2p ^- ( 2y'A - yA^- 

it V -T 



p-i 



x 



w—l 

E 

i=o 



(ID 



^5 (p)P (n-l-2p)i2«- 2 -^Aj] 



p=2 



(p-i)i(-ir 

' — 1 — w)! 



- 1 - j)! ^A 2 / 2(w- j) + 1 



'• 7 



x ^ F 



p=2 
w—l 

+ E 



10 



-J w\(p — 1 — iy)! 
2jR' 2j - x A 1 -^ 



A ' ^ j!(w-l-j)! 2(w-j) + \ 
On the other hand, using Eqs. (I2.1ip and (14.411) . we obtain 



y 



2(w-j)+l 



SR. 



(A.23) 



[n/2] 



(n 2)X-2 g 5(p)(7l _ ! _ 2 p) J R«- 2 - 2 f(l - F) p i? /2/A 



2k? 



p=o 

n/2] 



+ ( !i -^ Eawf)fl „_ 1 _ 2p(l _ frl ^ 

1 p=l 



(A.24) 



Some useful cancellations allow us to derive 

[n/2] 

Hn-2 

R' 2i-;l 
p-i 



p=2 L ^ 

(P -!)!(-!)• 



+ E 



P-I 



^ iy!(p — 1 — w)\ 



A 



2 2/ F- i (^ + F)(-)+2F-W- 



A 



1?' 



(p - l)!(-l) w ^ (u; - 2jy 2 ( w -rt +1 (y 2 + 



+ E 



E 



^wl(p-l-w)!^ j!(w-l-j)! 



2(w-j) + l 



x 2(y 2 + F) 



7 )' + (W + ^)^} 



(A.25) 
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Expanding the first term, we finally obtain 

p AP A + P M M> 
R' 



v — 2 



-2p 



p-1 



X 



(p_ !)!(-!)■» 
^w\(p — 1 — w)\ 



w—1 



w 



3=0 



j\(w -l-j) 



2 (w-j ) + l 



By direct calculations, we can show 



and complete the proof. 



(A.26) 



(A.27) 



A. 4 Equation (14.541) in Lovelock gravity 



In this appendix, we show Eq. (I4.54p in Lovelock gravity. While we consider the spherically 
symmetric case (k = 1) in the maintext, we derive the equations for general k in this 
appendix. 



Let us start from the action in the form of flA.lOj) : 

2pk p - 1 R n - 2p - 1 y{N r A)' 



J M = — I d ya (P^ 



2 /SjH 

n p=0 



2pkP~ 1 N 
n — 2p 



( R n-2py> A -l + ( R n-2 P y( A -lV 



+ pR n - 2p - l \k p ~ 1 - (k-F) p - 1 j^(AN r y + A^NR')^ - Ay^J 

+ (n - 2p - 1) \ (k - Ff + 2pk p - 1 F\NAR n - 2 ~ 2p - 2pk p ~ 1 R n - 2p ~ 1 yA 



Using 



(A 



_ lv A(F' + 2yy') R" 



2R! 



R'A 



(A.28) 



(A.29) 
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for the second line and integration by parts for the first line, we obtain 

£ J d 2 ya {p) -2pk^k\N r + ^yiT-^yy 



(n-2)V^ ] 



2kI 



p=o 



+ piT-^jp- 1 - (k - Ff- 1 } { (N r + ^F> - Ay^ 
- pR n - 2p - 1 ^ : (k - Ff^F 1 + (n - 2p - 1) (k - Ff NAR n - 2 - 2p 



- 2pk p - 1 R n ~ 2p - 1 yA 



+(t.d.; 



(A.30) 



where we also used R' 2 A 1 = FA + y 2 A. Using 
NA 



R 



2K n 



p(k - F) p - X F' + (n-l-2p)(k- F} 



R 
(A.31) 



P M M 



(n-2)l£> ^ 
2«H 



»/2 A 
V- - V A nn-l- 

2_^a ip) —R 



2p 



p=0 



p(fc - - (n - 1 - 2p)(k - Ff 



(A.32) 



which can be obtained from Eqs. ( 12. lip and ( I4.4ip . we obtain 



*-M 



NA 



d 2 y[P M M+—M' 



2k1 



(k) [n/2] 



jZjd 2 m P ) -2pk p - l A^N r + ^yR- 2 ^y)> 



+ pE 1 - 2 ?- 1 ^- 1 - (k — F)P- 1 j(^N r + ^F' - 2pk p - 1 R n - 2p - 1 yA 

+ (t.d.). (A.33) 
Using Eq. ( 1A.2j) for yAF/F together with integration by parts and R/R = Ny/R' + N r , 



pk^R^^Ayj + jr(k- F) p d t {R n ~ l - 2p ) 
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we rewrite the above expression as 



NA 



d 2 y(P M M+—M' 



(n - 2)K 



(As) [»/2] 



n-2 



2< 



p=o ^ 



-2pF- 1 A^(^- 2p - 1 2 /) / 



+ piT" 2 ^ 1 jp- 1 - (Jfe - F) p - 1 |^^F / + ^(ife - Ffd t {R n - 1 - 



2p\ 



+ pk p - 1 \2yAd t (R n - 2p ' 1 ) + — 

[ n — 2p 

Replacing the logarithmic term by 

R! + 



d t (R n - 2p )d x In 



R' + yA 



d x \ In 



R'-yA 



R' - yA 
-2R"yA + 2(yA)'R' 



■(t.d.). (A.34) 



FR' 



FA 2 

-F' + — F 

1/ 



we finally obtain 



(A.35) 



I M - I d 2 y(P M M 



2 k 2 



p=o 



(n - 2p - l)(Jfe - F)F' - pRF' 



R'F 



+(t.d.; 



/^M' + (,d. 



(A.36) 



and completes the derivation. 



B Boundary condition at spacelike infinity 

In this appendix, we confirm that under the transformation from {A, Pa; R, Pr} to {M, Pm\ S, Ps} 
the boundary term and total variation in the (14. 48ft are finite with suitable asymptotic fall- 
off rates of the ADM variables consistent with asymptotic flatness. 

We consider the following behavior of the ADM variables near spacelike infinity: 

A~ +0(aT ei ), (B.l) 

N r ^ N™(t)x- {n - 3)/2 ~ e2 , (B.2) 

A~ 1 + A 1 (t)x- {n - 3) - t3 , (B.3) 

R~ x + Fi(t)x- ( "- 4) - £4 , (B.4) 
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where e\ is positive and €2-64 are non-negative numbers, and require the following three; 

(I) the canonical transformation from {A, Pa; R, Pr} to {M, Pm', S, Ps} is well-defined, 

(II) the Hamiltonian is finite in terms both of {A, P&; R, Pr} and {M, Pm', S, Ps}, and 

(III) the Misner-Sharp mass is non-zero finite at spacelike infinity M ~ M°°(t). 

We will see below that these requirements are fulfilled for e\ > 0, e 2 > max[0, — (n — 5)/2], 
e 3 = 0, and e 4 > max[0, — (n — 5)], which are adopted as the boundary condition in the 
present paper. 

First we see the following integrated Liouville form (14.48j) : 

/oo poo poo 

dx(P A 5A + P R 5R) - / dx(P M 5M + P S 5S) = 5 r]dx + [( (B.5) 
00 J —00 J —00 

where P A , Pr, S, P s , M, P m , V, and ( are defined by Eqs. (@2D, (BSD, fTCTj) . flH| . 
( 14.42)) . ( 14.43)) . (14. 49 p . and ( I4.50p . respectively. For the well-definedness of the canonical 
transformation, two conditions must hold at spacelike infinity; (i) ( vanishes, and (ii) the 
integrands in the left-hand side and rj converge to zero faster than 0{x~ 1 ). The second 
requirement ensures the finiteness of the integrals. 

Let us see the asymptotic behavior of the momentum conjugates. Near spacelike infinity, 
Pa and Pr behave as 

P A ~ _ {n ~^ An - 2 N^ (r iX 1 ^ - N^-^-A , (B.6) 

Pr * - {n ~^ An - 2 N-' ^(t) (-^ + e 2 y n - 5 ^ + A^x^ + (n - 3)R lX ^ , 

(B.7) 

with which we obtain 

P A 5A ~0{x- (n - 4) - e3 - ei ) + 0{x- {n - 3)/2 - e2 - £3 ), (B.8) 
P R 6R ~0(x~ (n - 3)/2 - €2 - e *) + 0(x~ {n - 4) - e *- e4 ) + 0(x- (n - 4) - 2e4 ). (B.9) 

Hence, e 4 > and e 2 + £3 > are required for n = 5. For n = 4, the requirement is 
e 4 > 1/2, e 3 + e 4 > 1, e 2 + e 4 > 1/2, and e 2 + e 3 > 1/2. 

For the next check, we see the behavior of P, defined by Eq. (13.101) . Using 

y ~N^{R lX - {n - A) - ei - N™x- {n - 3)/2 - e2 ), (B.10) 

where y is defined by Eq. (13.121) . we obtain 

P -N-ZRiX-^i-Rtx-^- 264 + 2N™x- {n - 3)/2 - e2 - ei ) 

+ 1 - 2A 1 aT ( ™- 3) ~ e3 - N^° 2 x~^ n ~ 3 ^~ 2e2 - 2(n - 4 + e 4 )Pix" (n " 3) " e4 , (B.ll) 
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which converges to 1. For the condition (III) required, F must behave near infinity as 

F~ 1 -Fi(t)aT (n-3) (B.12) 
and then, for M ~ M°°(t), Fi is identified as 

m ^ 7 ^a. ( b,3) 

(n - 2)A n -2 

Since we have already required (n — 4) + 2e 4 > 1 and (n — 3)/2 + e 2 + e 4 > 1 in the previous 
argument, the first line in (IB.llj) converges to zero faster than 0(x~^ n ~ 3 ^). Therefore the 
condition (III) requires e 2 e 3 e 4 = 0, where Fi(t) is determined depending on the cases. 

Using Eq. (IB. 10p and M ~ M°°(t), we obtain the asymptotic behavior of PmSM as 

P M SM ~ -7V- 1 5M 00 ( J Ri^- (n - 4) - e4 - A^ r °°a;- (n - 3)/2 - e2 ). (B.14) 

This provides additional requirements; e 2 > for n = 5 and e 4 > 1 and e 2 > 1/2 for n = 4. 
Note however that if one uses the momentum variable Pm appropriate for asymptotically 
PG slices, then there is an extra term on the right hand side of ( 1B.14j) that cancels the 
second term. This in turn allows the choice e 2 = for all spacetime dimensions as required 
by the PG slicing. 

The conditions for e 2 and e 4 obtained up to here are summarized as e 2 > max[0, — (n — 5)/2] 
and e 4 > max[0, — {n — 5)]. We will see below that the requirements (I) and (II) are fulfilled 
under these conditions. 

Using the following asymptotic expansion of Ps] 

P s ~ - ( n ~ 2 J An - 2 N^ ^A lX - e3 + (n-4 + e 4 ) J R 1 x- e4 ^+C(a;- (n - 1)/2 - e2 ^) ) (B.15) 

where €m is some positive number defined by the next leading-order of M as M ~ M°° + 
0(x~ eM ), we obtain 

P S 5S ~ - 5R 1 ^—^=^ N- 1 (k lX -^-^ + (n - 4 + e 4 )R lX -^~ 2 ^ 

+ 0(x~ 3 ^ 3)/2 - e2 - e4 - eM ). (B.16) 

Under the present conditions this converges to zero faster than 0(x~ 1 ). Here too one 
should keep in mind that with the choice e 2 = needed for PG slicings, there is an extra 
total variation arising from the transformation from Pm to Pm- This extra term cancels the 
second term on the right hand side of ( IB. 161) rendering i] finite for e 2 = in all dimensions. 

Next let us evaluate ( and 77. We write the logarithmic term in the expressions (I4.49P and 
(1450]) as 



R'-yAJ \l + W 
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where 



W :-- 



{n-2)A n - 2 R n - 3 R r 
Using the fact that W converges to zero as 



W ~ -A^ 1 ^i? ia ;- (n - 4) - e4 - N^x' {n - 3)/2 - e2 ^j . 



we can evaluate the logarithmic term as 

\R'-yAj 3 
Using this, we can show that ( converges to zero under the present conditions as 



C ~0(aT (n - 5) ~ 2e4 ) + 
On the other hand, r\ is evaluated as 

(n - 2)A n - 2 R n ~' i R' 



0(x 



-{n-5)/2-e 2 -ei\ 



(B.18) 



(B.19) 



(B.20) 



(B.21) 



V 



2k^ 
[n - 2)A n - 2 x 
3^! 



2lU + ln 



1 - W 
1 + W 



n— 3 



-w 6 



^ 0(a; -2(n-4) + l-3e 4) Qr ^ n -3)/2-^ 

and converges to zero faster than 0(x~ v ) under the present conditions. 



(B.22) 



Lastly, let us check the well-definedness of the Hamiltonian dx(NH + N r H r ), where 
H r and H are defined by (14.71) and (14. 8p . respectively. It requires that the fall-off rates of 
NH and N r H r are faster than 0{x~ 1 ). From the following asymptotic expressions; 



'1 - eAR lX -^ - 



n — 3 



^2 



(B.23) 



R'P, 



R 



(n-2)A>-2 N i / 1 _ , _ 4 + )a .-(»--»)-. , 



n — 3 



+ e 2 j x (n " 5)/2 ~ e2 + Ai(*)af 63 + (n - 3)^ 



.r 



-«4 



(B.24) 



we see that the dangerous terms of the order x*™ 5 - ,//2 62 in if r = —AP^ + R'Pr are canceled 
out and obtain 



N r H r ~0{x 



-(n-3)/2-e 2 -e 3 \ 



+ 0(x 



-(n-3)/2-e 2 -e 4 \ 



(B.25) 
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This is faster than 0{x 1 ) under the present conditions. Next let us see the behavior of 
NH. Using the followings asymptotic expansions; 

ad ^ {n - 2)A n - 2 _ N _ x 



(l + A ia ;- {n - 3) - e ^ ^x 1 - 64 - N™x {n - 3)/2 - €2 j , (B.26) 

RPr * - {n ~ y** 1 - 2 N- 1 (x + fl ia r<"-*)-*) 



B.27) 



we obtain 

Pa f "-3 



i? 



n-2 



+ C(x- (n - 4) - 2e4 ). (B.28) 



This shows that the first line of the following super-Hamiltonian; 



H _ J RPr -^1a Pa 



"n^ A 

{n-2)A n - 2 R n - 2 V^ " 2 



(n 2)A n - 2 / _ jR n-3 (jR / A -ly + !^ AjR n-4 (1 _ , ( y^f) ) 



converges to zero faster than 0{x~ 1 ) under the present conditions. On the other hand, the 
second line is evaluated as 

-iT-^P'A- 1 )' + ?l^Air- 4 (l - A- 2 P 2 ) ~ Oix- 1 - 63 ) + C^aT 1 " 6 *), (B.30) 

which also converges faster than 0(x~ l ). As a result, NH converge to zero faster than 
0{x~ 1 ) and hence the Hamiltonian is well-defined. 

We also check the well-definedness of the Hamiltonian with a new set of variables: dx(N M M'+ 
N S P S ), where N M and N s are defined by (14.561) and 04.571) . respectively. Using Eq. (IB.lOp . 
we obtain A^ A/ ~ N^. Combined this with M' ~ 0(x~ 1 ~ £M ), it is shown that N M M' con- 
verges to zero faster than 0{x~ x ) under the present conditions. On the other hand, using 
N s ~ C(x- (n - 4) - £4 ) and Eq. flB~15l) . we obtain 

N S P S ~C(x~ (n - 4) ~ e3 - £4 ) + C(x" (n - 4) - 2e4 ) + 0(x- 3{n - 3)/2 - €2 ~ ei ~ eM ). (B.31) 

This also converges to zero faster than 0(x~ l ) under the present conditions. 
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